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Dispersion interactions such as the van der Waals interaction between atoms or molecules derive 
from quantum fluctuations of the electromagnetic field and can be understood as the exchange of 
virtual photons between the interacting partners. Any modification of the environment in which 
those photons propagate will thus invariably lead to an alteration of the van der Waals interac¬ 
tion. Here we show how the two-body dispersion interaction inside a cylindrical waveguide can be 
made to decay asymptotically exponentially, and how this effect sensitively depends on the material 
properties and the length scales of the problem, eventually leading to the possibility of controllable 
interactions. Further, we discuss the possibility to detect the retarded van der Waals interaction by 
resonant enhancement of the interaction between Rydberg atoms in the light of long-range potentials 
due to guided modes. 

PACS numbers: 34.20.Cf, 34.35.+a, 42.50.Ct, 42.82.Et 


I. INTRODUCTION 

Dispersion interactions such as the van der Waals 
(vdW) interaction are best known for generating attrac¬ 
tive potentials between electrically or magnetically neu¬ 
tral yet polarizable quantum emitters as a direct con¬ 
sequence of the quantized nature of the electromagnetic 
field. While seminal work studied either such interac¬ 
tions between two quantum emitters embedded in a bulk 
medium m or the (Casimir-Polder) dispersion inter¬ 
action of a single quantum emitter with a macroscopic 
body mil], a question of interest is how two-body dis¬ 
persion potentials can be modified by the presence of 
macroscopic external boundaries that affect the mode 
structure of the electromagnetic fields and, hence, also 
the vacuum fluctuations mu. Very strong effects can 
be expected in waveguides where well-dehned modes in¬ 
duce long-range correlations. For a single emitter in a 
cylindrical waveguide near a cavity resonance, Ref. m 
predicted strong enhancement of the single-particle-wall 
interaction at certain resonant radii. Recently, the two- 
body dispersion interaction in an idealized metallic rect¬ 
angular waveguide was shown to decay asymptotically 
exponentially with distance m if a lower mode cut-off is 
present, whereas the fundamental mode of a transmission 
line provides a signihcant enhancement of the dispersion 
potential [T^ . 

Here, we consider the case of hollow cylindrical tubes 
(capillaries) of a metallic or semiconducting material, or 
dielectric cylinder waveguides (optical fibers) that con¬ 
tain a pair of quantum emitters, as depicted in Fig.^ 
and apply the approach of macroscopic QED to gain a 
deeper understanding of the underlying mechanisms and 
study the impact of material properties beyond the ap- 
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proximation of ideal conductors. We present numerical 
results that agree well with analytic calculations. Our 
results show how a modification of the boundary con¬ 
ditions imposed by a realistic material can be used to 
modify the van der Waals potential significantly. The 
drastic difference between dielectric and metallic bound¬ 
ary conditions make the use of phase-change materials 
particularly interesting. 

We further discuss the role of resonant potential con¬ 
tributions which arise from the emission of real photons 
by initially excited atoms. In a waveguide environment, 
these may bring a discussion of the still elusive far-held 
potential into reach. Owing to the rapid decrease of the 
vdW interaction with the emitter separation, detection 
of the retarded potential is generally exceedingly difh- 
cult. A possible alternative to use highly excited Ryd¬ 
berg states with their vastly enhanced vdW interaction 
mnii will be discussed. A detailed review of the dyadic 
Green function in a waveguide environment is given in the 
Appendix. 


II. WAVEGUIDE-ASSISTED DISPERSION 
POTENTIAL 

A. General approach 

The vdW dispersion potential is obtained from the mu¬ 
tually induced dipole-dipole interaction within fourth- 
order perturbation theory. We obtain the general ex¬ 
pressions within the framework of macroscopic QED 
n na HD, where the mode structure of an arbitrary 
dispersive and dissipative environment is encoded in the 
classical dyadic Green function Q{r,r',uj). 

Evaluating the dispersion potential near an interface, 
one finds single-body energy shifts which depend only on 
the position of a single emitter laii]. Here, we focus on 
the irreducible two-body contribution that depends on 
the positions of both emitters. This dispersion potential 
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FIG. 1: Sketch of the system: two emitters are placed on the 
axis of a cylindrical waveguide of radius R, separated by a 
distance a. The core medium £2 is embedded in a low-index 
(or metallic) medium £1 to form a dielectric (or conducting) 
waveguide. 


is given by the well-known expression [ilia [III nil 
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For two-level emitters with resonance frequencies fin and 
transition dipole moments in the ground state, we 
have the polarizability tensor 
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We consider two identical particles with a preferred po¬ 
larization axis. Denoting by a and Q the projections of 
the polarizability and the Green tensor onto this direc¬ 
tion, we obtain scalar expressions that can be expressed 
at imaginary frequencies as 

u = -hf ^-^a^ii^)g^{ri,r2/i0 ■ (3) 

Jo 27 r 

We now turn to the setup sketched in Fig.j^ where two 
identical particles with either purely axial or radial polar¬ 
izability are placed on the axis of a cylinder of radius R 
and separated by a distance z. In the presence of bound¬ 
aries, the scattering ansatz is employed Q = 
where and refer to the bulk and scattering parts 
of the Green tensor, respectively. Explicit expressions for 
the Green tensor in a cylindrical waveguide are compiled 
in Appendix [^ This leads to the decomposition of the 
vdW potential into three terms Bng, 


g2 ^ (g(O)y ^ ( 4 ) 

-f -h . (5) 


Each of these three terms has a clear physical meaning: 
the first ) describes a loop of freely propagating pho¬ 
tons that dominates in the near-held regime of distances 
z X, R, where it recovers the two-body potential in a 
homogeneous medium, i.e. the usual dipole-dipole near- 
held coupling with a distance scaling. The second 


FIG. 2: Contributions to waveguide-mediated vdW cou¬ 

pling between two emitters inside a subwavelength hollow 
metallic cylinder (black curves) or dielectric cylinder waveg¬ 
uide (undoped silicon in vacuum, green curves) with radius 
R = 0.8c/D. a) Radial dipole orientation and b) axial dipole 
orientation. The arrows indicate the onset of retardation at 
z-^fe^Sljc Ri 1. 


term (C/^®'^)) is due to loops formed by two scattered pho¬ 
tons, and depends on the boundary conditions. At short 
distance, it approaches a constant value. The third term 
(?7(cross)) consists of two-photon loops that involve only a 
single scattering event, and has a subleading scaling 
at short distance. 

In the retarded regime z ^ A, each of the three terms 
scales as r“^. This is clearly seen in Fig.[^) and b), 
where the three contributions are plotted for perfect 
metallic boundaries and for a dielectric waveguide. The 
two panels correspond to radial and axial dipole orienta¬ 
tion, respectively. The material properties of the cylinder 
walls do not qualitatively change the scaling of the indi¬ 
vidual contributions to the dispersion force potential. In 
fact, a change in the boundary condition will mainly af¬ 
fect the peculiar properties of the discrete guided modes 
that can propagate in the waveguide. For all conduc¬ 
tors, material properties have relatively little impact on 
the retarded potential which is dominated by the fluc¬ 
tuations of evanescent modes. We can expect that, in a 
scenario where one of the emitters is initially excited, the 
exchange of real photons may give rise to resonant con¬ 
tributions to the van der Waals interaction resulting in a 
much stronger dependence on the boundary conditions. 
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B. Perfectly reflecting cylinder 


We first calculate the interaction potential in a per¬ 
fectly reflecting waveguide using the scattering decom¬ 
position. In the far field (z ^ A, i?), the scattered Green 
tensor in an ideally reflecting waveguide exactly ap¬ 
proaches the homogeneous one except for a sign change, 
which is mirrored in the contributions to the dispersion 
potential shown in Fig.[^) and b) for radial and ax¬ 
ial dipole orientation, respectively We thus find from 
Eq. ([I I that the cross-term nearly cancels the other two 
contributions, except for a small remainder with expo¬ 
nentially suppressed distance behavior, visible in Figs.[^ 
a) and b) that show the potential normalized to the value 
in free space for different cylinder radii. This surprisingly 
exact balance lies at the origin of the exponential behav¬ 
ior of the retarded potential, similarly encountered in the 
case of rectangular waveguides m 

In order to see the exponential decay directly, it is help¬ 
ful to use an alternative approach to the dyadic Green 
function that is better suited for perfectly conducting 
cylinders. As a consequence of Maxwell’s equations, the 
tangential components of electric field have to vanish on 
the cylinder surface. Hence, the vector wave functions 
used in the expansion of the Green function are [18] 


Men,ih)=Wx 
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with fi = qnm/R and A = Pnm/R- Here, the numbers 
Pnm denote the mth root of the Bessel function of order 
n, Jn{Pnm) = Oj and qnm the mth root of the derivative 
of the Bessel function of order n, J^iqum) = 0. 
Expanding the Green function gives for 

g{r,r ,w) =-r^A(r-r ) 
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FIG. 3: Dispersion potential for two emitters in a perfectly 
reflecting hollow cylinder waveguide normalized to the free 
space potential (dotted line) with axial (a) and radial (b) 
polarizability for cylinder radii RQ/c = 0.4, 0.8,1.6 (corre¬ 
sponding to 50nm, lOOnm, 200nm at A = 780nm). The dots 
indicate the asymptotic far field behavior of Eq. (181 to the 
axial component, and exponential fits to the radial compo¬ 
nent, respectively. 


On the cylinder axis and setting z > z' = 0, only very 
few of the terms survive, and the result can be written 
as 
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The van der Waals interaction between two dipoles in 
a perfectly reflecting cylinder possesses only two intrin¬ 
sic length scales, the radius R of the tube and the dipole 
separation z. The limit z/R —>■ 0 corresponds to an in¬ 
creasingly large cylinder radius R or, alternatively, very 
close dipoles. In both cases, the van der Waals interac¬ 
tion in free space is recovered as follows. In the limit of 
large arguments, the Bessel functions become 



(13) 


and their zeros take the form pnm — {jn + f ~ i) ■ 
Hence, the difference between adjacent zeros is Ap^m = 
Pn,m+1 - Pnm — TT and thuS AA = Apnm/R TT/ R. In 
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this way, we find for example 
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which exactly reproduces the (n = 0)-contribution to the 
free-space Green function [la¬ 
in the opposite limit, zjR ^ oo, the cylinder radius 
is small compared to the dipole separation. In this case, 
we can approximate the exponents as 

k\z = x/{kRf - pI^ ~ , (15) 


such that for an axial dipole the dominant contribution 
to the Green function reads 


Gzz{r,r',uj) ~ 


Poi 


27rfc2i?3 J^(poi 
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where only the lowest (TMoi) mode with poi — 2.405 
contributes. The van der Waals potential derived from 
this Green function is thus 


Uzz{z) 
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which is seen to decay exponentially with increasing 
dipole separation. Gompared to the van der Waals in¬ 
teraction in free space in the nonretarded limit, the ex¬ 
ponential suppression is 


Uzzjz) 

ui^’iz) 


‘^Poi 

JfiPOl) 




(18) 


Good agreement of the scaling is seen in Fig.[^). The 
change of sign encountered in this component leads to 
an enhancement of the potential U^z for axially polar¬ 
ized particles above the limit in a homogeneous medium 
at short distances. In contrast, the potential for ra¬ 
dial polarization is always below the free-space value 
U±(z) < U^\z). For this component the logarithmic 
plot suggests a different far-field scaling U±{z) oc e~^z~'^ 
[fits in Fig.j^D)]. 


C. Conducting cylinder 

We now assess numerically the impact of realistic ma¬ 
terials. A concise representation of the Green function 
for a cylinder based on the scattering decomposition is 
given in App. The general boundary value problem is 
expressed in terms of reflection matrices that contain the 
dielectric functions, e.g. for a Drude metal and vacuum 

ei(w) = l-, £2 = 1. (19) 

a;(a; -I- 17 ) 
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FIG. 4: Dispersion potential for two optical emitters in a 
hollow conducting cylinder {fl — 27rc/780nm, Rc/fl = 0.8) 
with axial (a) and radial (b) polarizability, normalized to the 
potential in free space. The curves correspond to a hollow 
waveguide with perfectly reflecting boundaries (black) or the 
material boundaries for gold (gray), VO 2 (blue), and doped 
silicon. The red dots indicate the asymptote [Eq. @1 to the 
axial potential and suggest a weak dependence on material 
parameters. The red squares and diamonds are obtained from 
an exponential fit to the radial potential with a significant 
material dependence. 


For gold, the relevant parameters are hiOp = 8.5eV, 7 = 
5 X IQ-^Wp [201. 

At optical wavelengths (A = 2^01^1 = 780nm, i? = 
0.8c/fl = lOOnm), we find that the impact of imperfect 
boundary conditions on the radially oriented dipoles is 
much stronger than on the axial one, see Fig.[^ Perfect 
boundaries and good conductors agree closely. The ex¬ 
ponential fit to the potential for radially aligned dipoles 
agrees with an effective radius i?off = ISOnm, while the 
axial one is well-described with the asymptote of Eq. (18 1 . 

In order to understand the impact of a limited con¬ 
ductivity, we compare to the results obtained for hollow 
core fibers made out of weakly doped silicon (impurity 
concentration 1.3 x 10^®cm“^) and metallic vanadium 
dioxide (VO 2 at T > 340K), which can be described 
by Lorentz-Drude models with the parameters given in 
Refs. [HJ [22]. A numerical evaluation of the distance- 
dependent potentials (normalized to the bulk potentials) 
is shown in Fig.|^ As the conductivity is reduced, the 
boundaries become more penetrable, and the cancella¬ 
tion of scattered and free contributions is no longer ex¬ 
act. We see that the potential for axially polarized emit- 
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FIG. 5: Normalized dispersion potentials for two opti¬ 
cal emitters inside a silicon cylinder in vacuum (fl = 
27rc/780nm, Rc/Q — 0.8) with axial (solid green curve) and 
radial (dashed green curve) polarizability. The results for a 
perfectly conducing hollow tube (black curves) are given for 
comparison. The normalization refers to the bulk potential 
in the core medium. 

ters inside a conducting V02-capillary begins to deviate 
slightly from the perfect conductor limit, whereas for the 
doped silicon capillary, the far held potential recovers the 
z“^-behavior. Again, the impact is much more drastic for 
the radial component. 

An interesting question is whether the superconduct¬ 
ing phase transition can lead to a signihcant modihca- 
tion of the dispersion potential. For an impure niobium 
sample the parameters in Eq. ( |19[ ) should be chosen as 
fuMp = lOeV, 7 = 6.5 x |23j . A minimal - yet rea¬ 

sonable [24] - Meifiner-London model for the supercon¬ 
ducting state far from criticality is obtained by setting 
7 —>■ 0. Here, the Meifiner length A^ = c/wp « 20nm 
describes the penetration of helds into the surface. At 
optical frequencies, the superconducting phase transition 
does not have a great effect, as the penetration into nor¬ 
mally conducting niobium, described by the skin depth 

(5 = y^ 27 c^/(WpW) ~ 5nm, does not differ greatly from 

the Meifiner-London length, so that both cases provide 
potentials similar to the one for gold in Fig.j^ For a 
lower transition frequency, e.g. for a Rydberg transition 
|n, 1 = n — 1) —>■ \n' = n + 1,1' = n) with n = 50 at 
H = 27r X 51GHz (A = 5.6mm), we find that the normal 
skin depth S « 500nm is much larger than the Meifiner- 
London length. Hence, the superconductor can be ex¬ 
pected to approximate even better an ideal reflector at 
GHz-frequencies. However, the variations are negligible 
with respect to the waveguide radius unless RQ/c 1. 

D. Solid dielectric cylinder 

Similar to the transmission-line structures considered 
in Ref. [12], a dielectric waveguide features a fundamen¬ 
tal mode without a frequency cut-off. Modifications of 
the results are expected due to the frequency-dependent 
mode profile in a dielectric waveguide. We assume solid 



FIG. 6: Normalized dispersion potential for two emitters in a 
silicon cylinder embedded in VO 2 (11 = 27rc/780nm, Rc/fl = 
0.8) with axial (solid cnrves) and radial (dashed curves) po¬ 
larizability. At T < Tc = 340K, the cladding is dielectric and 
provides a weak dielectric waveguide (green curves) with mod¬ 
erate enhancement Uzz{oo)/U^z { 00 ) « 2. Above the phase 
transition (T > Tc, black curves) the metallic cladding sup¬ 
presses the far-held potential. 


state emitters embedded in an undoped silicon cylin¬ 
der suspended in vacuum, and impose a Drude-Lorentz 
model [12] that agrees well with the results of a simple 
index medium (ei = 1, £2 = 13.5). 

The resulting dispersion potentials in this structure for 
the two polarizations are compared in Fig.[^ to the po¬ 
tential in a perfectly conducting hollow cylinder. For the 
radial dipole orientation, we observe an enhancement by 
a factor of 2 at intermediate distances zc/H « 1, similar 
to the effects found in the surface-assisted two-body po¬ 
tential [7] . This is connected with a change of sign in the 
cross term 7 /('=''oss) ^ ji^jicates the change from a bind¬ 
ing contribution at short distances to an antibinding con¬ 
tribution encountered at all distances for the conducting 
boundaries, see Fig.]^ For the axial dipole orientation, 
an enormous enhancement by three orders of magnitude 
with respect to the free space potential survives to large 
distances. 

In both cases, the power laws are recovered in the 
far field, indicating that the subtle cancellation at work 
for the metallic waveguides is perturbed, especially for 
the axial dipole components. This can be clearly seen 
from the single contributions in Fig.[^ The scattered 
contribution can be identified as the origin of the 

significant enhancement in the potential of axially polar¬ 
izable emitters. 


E. Metal-insulator transition 

The difference between a superconducting and a nor¬ 
mal metal were found to be small, as both impose sim¬ 
ilar metallic boundary conditions. However, the dras¬ 
tic qualitative difference between the potentials encoun¬ 
tered for dielectric waveguides and metallic ones, result¬ 
ing from the existence of a fundamental mode or a mode 
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cut-ofF, respectively, suggest that materials with a metal- 
insulator transition may have a significant impact. 

Vanadium dioxide (VO 2 ) is known to change from a di¬ 
electric to a metallic phase at Tc = 340K, and its impact 
on the Casimir interaction in coplanar cavities has been 
discussed before [22], based on Drude-Lorentz models 
of the dielectric function m- Interestingly, its refrac¬ 
tive index in the dielectric state is slightly lower than 
that of undoped silicon, so that a silicon core embed¬ 
ded in VO 2 may be switched from a dielectric guide to 
a metallic one by a small increase in the temperature. 
This affects both polarization components as is clearly 
visible in the normalized potentials in Fig.]^ Due to the 
low index contrast, we obtain a leaky waveguide and at 
T < Tc (green curves) the potential recovers closely the 
one in a homogenous medium. Still, the axial component 
is enhanced by roughly a factor of 2 . As T > Tc (black 
curves), the metallic boundary results in the exponential 
behavior already discussed in detail, so that the potential 
can be effectively switched off. Note that thermal correc¬ 
tions to the potential are negligible for optical transitions 
at the temperatures considered. 

We expect an improved performance in a coated fiber 
in vacuum, where a V 02 -layer in the insulating state 
would provide a situation close to a homogeneous di¬ 
electric cylinder with good waveguiding properties. This 
could recover the strong far-field enhancement encoun¬ 
tered before. To generate a sufficiently strong boundary 
condition in the metallic state, the layer thickness should 
however exceed the skin depth (d « 60nm at A = 780nm). 


III. OBSERVABLE FAR-FIELD POTENTIALS 


If one of the emitters is not in its fundamental elec¬ 
tronic state, the potential changes as a result of Forster- 
like processes HZmiMH] 

Ugeiz) = -U^dw{z) + Ures{z) . (20) 


The resonant potential arises from the exchange of real 
photons between the emitters via the guided mode. In 
metallic waveguides, this requires transition frequen¬ 
cies above the cut-ofF. In the single-mode regime, 
one can approximate the far field Green function by 
Q{z,uj) « i/3Im [0(0, D)] exp(inefffc| 2 |), where fi describes 
the emitter-wavegui de co upling efficiency and rieff the 
mode index (see App. A3). When bulk absorption is low, 
this results in interactions of infinite range, ultimately 
limited by the photon coherence 2 z « c/ 7 . 

Unfortunately, the form of the resonant potential is 
still disputed in the literature HSJ |25] and two forms are 
commonly obtained: 


=|d(i)|2a(2)(Ui)Re([0(ri,r2,Ui)]2) (21) 

OC COS^(nefffc|-2|) , 

a^"^(f^i) |0(ri,r2,Ui)|2 (22) 

OC const. 


This seems to originate in different assumptions about 
the system dynamics. In fact, the first result is com¬ 
monly encountered for a system initially prepared in a 
state |e) 0 |g) and left to free evolution, while the second 
arises from steady-state nonequilibrium configurations. 
Mathematically, the two forms result from different treat¬ 
ments of the double poles that arise in fourth-order per¬ 
turbation theory. The first result is obtained by treating 
the poles as principal values, the second makes use of 
Sokhotsky’s formula that also contains a d-function con¬ 
tribution. On the basis of perturbation theory alone, the 
correctness of one or the other form of the resonant po¬ 
tential cannot be decided, as perturbation theory makes 
no statement as to how the poles should be circumvented. 
While the first case would yield an oscillating potential 
and force - as in free space -, the second potential is con¬ 
stant and, in consequence, force free. Arguments relating 
the resonant vdW potential between two dipoles to the 
dilute limit of the Casimir-Polder interaction between a 
dipole and a macroscopic body (known also nonperturba- 
tively [29]) suggest that correctly describes systems 
left to free evolution [5D]. A dynamical approach that 
involves solving the coupled atom-field dynamics within 
the framework of macroscopic QED supports this view 

m- 


More importantly, however, the nondecaying distance 
dependence may allow for the observation of the far-field 
dispersion potential between emitters, which has so far 
remained elusive due to the quickly decaying behav¬ 
ior in free space. At this point one might be tempted 
to try and use highly excited Rydberg atoms with their 
exaggerated vdW interaction to detect the retarded in¬ 
teraction. However, a straightforward inspection of the 
relevant scaling laws with the principal quantum number 
n [ 22 ] shows that this is not advantageous: the dipole mo¬ 
ment matrix elements of neighboring dipole-coupled Ry¬ 
dberg states scale with [dj oc n? and the associated tran¬ 
sition frequencies as D ex n~^. Hence, the polarizability 
is found to scale as a (x . As the boundary between 
nonretarded and retarded regimes is set by the transition 
frequency as z ~ c/0, we immediately find that the vdW 
potential U^dwiz^) at this distance is in fact suppressed 
by a factor oc compared to the ground-state poten¬ 
tial. Similarly, resonant Forster processes that behave as 
Uces{z) = C 3 Z~^ with C 3 oc are suppressed by a factor 
oc This is due to the fact that the benefit of a larger 
interaction strength is strongly outweighed by the rapid 
increase of the length scale at which the retarded regime 
sets in. 


However, the resonant potential that requires prop¬ 
agating photons that are guided along the waveguide 
structure does not have the limitations as in free space. 
In this case one would expect to reach the retarded limit 
between two Rydberg atoms with an appreciable inter¬ 
action strength. 
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IV. CONCLUSION 


homogeneous medium Si, and decompose 


We have studied the vdW dispersion potential between 
two emitters placed inside a realistic metallic or dielec¬ 
tric waveguide with a subwavelength diameter. An expo¬ 
nentially decreased interaction potential was found for a 
perfectly reflecting boundaries, similar to the results for 
rectangular cross sections m- We found that the range 
of the potential can be modified by tuning the conduc¬ 
tivity of the capillary wall. Fluctuations mediated by the 
fundamental mode in a dielectric waveguide, in contrast, 
can strongly enhance the far-held potential. This was 
previously unsettled due to the 3D nature of the guided 
modes [H]. The drastic qualitative difference between 
the potentials could allow to realize switchable disper¬ 
sion potentials between solid-state emitters embedded in 
a dielectric waveguide core coated by a cladding featur¬ 
ing a metal-insulator phase transition. Finally, resonant 
contributions to the potential, arising from real-photon 
exchange via the guided mode, could bring the retarded 
dispersion interactions into the reach of experiments. 
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• ^ 77 

/ ^(2-Mx (A4) 


rMMMe^{h, 772) 0 772) 

-b rMNNe^{h, 772) (g) 7 V'e„(-/l, 772) 
± rMNMo^{h, 772) g) N'e^{-h, 772) 

± rNMNo^{h, 772) g) M'e^{-h, 772) 


Summation over upper and lower indices and global signs 

is assumed, e.g. ±Ai(^B'k = +Ai®Bj — Ak®Bi. The 
3 I 

cylindrical vector wave functions are 


We gratefully acknowledge fruitful discussions with 
Francesco Intravaia, Simen A. Ellingsen, Xuewen Chen, 
loannis Chremmos, Helge Dobbertin, and Vahid San- 
doghdar. This work was partially supported by the DFG 
(grant no. SCHE 612/2-1). 


Appendix A: Dyadic Green tensor in a cylinder 


1. Bulk Green tensor 


Me {h,r]) = V X 


Jnivp) (ncj)) 






Jn{r]p) {n(l)) 


(A5) 


(A6) 


where upper and lower cases refer to even and odd func¬ 
tions in (j). The function Jn{x) denotes the Bessel func¬ 
tion of the hrst kind, and a prime refers to the second 
spatial argument. 

In the case of a single interface, we generalize the forms 
given in Ref. [M], to obtain the concise expressions 


The Green tensor describes the field generated by an 
oscillating dipole source. In a bulk environment the held 
generated at a distance along the 2 :-axis can be expressed 
in the compact form m, 

g^°\r,r',uj) = -^^{l-ikz + k^z^) (AI) 

0(°)(r,r',w) = ^|p^(2-2ifcz) . (A2) 

Here, k = elzP' jt? denotes the propagation constant in 
a medium with permittivity e and = Qyy^ = G±'^- 
The form is useful for efficient numerical treatment. It 
can also be expressed in terms of cylindrical vector wave 
functions [3S]. 


2. Scattered Green tensor 

We consider a nonmagnetic dielectric cylinder of ra¬ 
dius R (dielectric function £ 2 ), surrounded by another 


r-(") — _ 


H, 


( 1 ) 


{x) A- 

X 




“ Jnix) A + C 
A = -n^(Ruj/cf(hRf(£i - £2 f, 


(A7) 
(AS) 

C = xfx2[h(xi)x2 - j(x2)xi][£ih(xi)x2 - £2j(x2)Xl], 

(A9) 

BmM = XiX2[h(Xi)x2 - h(x2)xi][£ih(xi)X2 - £2j(x2)xi], 

(AIO) 


BjvN = XiX2[h(xi)x2 - j(x2')Xi][£ih(xi')X2 - £2M^2)xi], 

(All) 

Bmn = -A + mx^X2(Rk2)(hR)(£2 - £i)lh(x2) -j(x2)], 

(A12) 


where Xi = R^/k'f — h'^. We abbreviated h{x) = 
^\nH^\x), j(x) = ^InJ("')(a;), where h['^\x) de¬ 
notes the Hankel function of the hrst kind. It may be 
useful to notice that the diagonal rehection matrices are 
even functions of h, while the off-diagonal ones are odd. 
By reciprocity, xmn = xmm, and causality requires that 
g{r,r’,-u3*) = g*{r,r\uj). 














For ideally reflecting surfaces (ei —>■ — cxd) the reflection 
matrices take the simple form [34] 


(n) _ („) _ H^\x)h{x) 

Ux) ’ Ux) j{x)’ 

^MN = = 0 • 


(A13) 

(AM) 


On the cylinder axis {p,p' 0) only the lowest 

multipole waves with n — 0,1 contribute. We choose 
(j) = (f)' = 0 and z' = 0 for simplicity, so that the Green 
tensor depends only on the axial displacement z. This 
gives 
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^1.2' NN 
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(A15) 


p2^NN 

"-2 


The integral can be immediately performed at imagi¬ 
nary frequencies. At real frequencies, parity arguments 
allow to map the integrand to w G [0,oo]. A numerical 
integration can be performed along a path in the com¬ 
plex h-plane that avoids the region h G [.y^w/c, 


from below and is slightly shifted into the upper plane 
elsewhere for improved convergence. 


3. Guided mode contributions 


The fundamental guided mode corresponds to the pole 
of the reflection matrices at a value h = neffLojc, which 
can easily be found numerically. The residue theorem 
is applied to evaluate the pole contribution. This cor¬ 
responds to replacing ppf —>■ i 7 ri 5 (h — h) gf^+P/dh ™ 
Eqs. (A15) and (AI 6 I. The pole contribution provides 
dominates the far field, so that 


Q{z,uj) « i/3Im [^( 0 , 0 ;)] exp(meff|- 2 |<j-'/c) (A17) 

and has a harmonic spatial dependence |35j . In lossy en¬ 
vironments, where the poles are removed from the real 
axis, their contribution can be found by a numerical con¬ 
tour integration or a Lorentzian fit to the integrand [36j . 
Note that the group velocity of homogeneous propaga¬ 
tion due to differs from the one of the guided mode 
and may result in a small beating of the amplitude at 
short and intermediate distances. 
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